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C^ ■ Explicit formulae are given for the saddle connection of an integrable family of 

standard maps studied by Y. Suris. When the map is perturbed this connection is 

destroyed, and we use a discrete version of Melnikov's method to give an explicit 

formula for the first order approximation of the area of the lobes of the resultant 

turnstile. These results are compared with computations of the lobe area. 
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1 Introduction 

Standard maps are area-preserving diffeomorphisms of T x R given by 

ie', r') = fie, r) = ie + r + v'ie),r + v'ie)) , (i) 



*A version of this paper was originally written in 1996, but never published. Subsequently, the methods 
of this paper have been used and extended, for example by Delshams and Ramirez-Ros [y . 
UDM was supported in part by NSF grant number DMS-9971760. 



where the potential, V, is periodic, V{6 + 1) = V{6). In this paper we study the standard 
map fs given by (|l|) with the potential^ 



Vs{9) = -- [ dt tan'^ (- 
^ Jo VI 



(5sin(27rt) 



(2) 



+ 6 cos (27rt) 

Suris [0, ^ showed that fs is integrable with integral 

Is{d, r) = cos nr + 6 cos tt(26 — r) , 

i.e.. Is o f^ = Is. Contours of Is are shown in Fig. |l]. The map is integrable for any 6; 
however, we will consider the case < (5 < 1, as the topology of the saddle connections 
changes at 5 = 1. For < 6 < 1 the map fs has hyperbolic fixed points at ^^ = (— |, 0) 
and Zb = (|,0), that are connected by two saddle connections, forming the upper and 
lower separatrices of the fixed point resonance. 




Figure 1: Some contours of Is, for 6 = 1/3. 

Standard maps are examples of monotone twist maps (for a review, see P]) defined 
by the geometrical twist condition that vertical lines tilt to the right upon iteration, 
^ > 0. Such maps have Lagrangian generating functions, S{9,9'), which generate the 
map implicitly through the equation dS{6, 9') = r'dO' — rdO, or equivalently 

r = -diSie,9'), r' = d2S{6,9'), (3) 

where the subscripts indicate derivatives with respect the first and second arguments of 
S. The twist condition implies that did2S < 0. We assume that our twist map has zero 
net flux, which is equivalent to S{6 + 1, ^' + 1) = S{6, 6'). For maps of the standard form, 
(0), the generating function is 



s{e,e') = -{e'~ef + v{e) . 

^ The potential can also be written 

V{9) = ^n [dilog(l + S)- dilog(l + Se^ 



(4) 



where the dilogarithm is defined by dilog(a;) = J^ °^_J dz. 



We will let Ss denote the generating function for the Suris map fs- 

Since the Suris map is integrable and has twist, Birkhoff's theorem [Q implies that 
the saddle connection between the equivalent points Za and Zb is the graph of a function 



,1 11 

2' 2J 

defines a circle diffeomorphism 



r = x{(^)^ for ^ e W = [— 2, 2]- The dynamics (|l|) restricted to the saddle connection 



e^h{e) = e + xi0) + v'{e). 

Our goal is to show that, when the Suris map is perturbed, the stable and unstable 
manifolds of the continuation of the fixed points intersect transversely. We use a version 
of Melnikov's method for twist maps [^ |^, ^, ||, ^ |1^, ll^]. The system is perturbed 
by adding to the original generating function any C^ periodic function P{6,6'). We will 
assume for simplicity that P does not move the fixed points. Then for small enough e, 
Ss + eP is a twist generating function, and the corresponding map fs,e has hyperbolic 
fixed points at Za and Zfy. As was shown in P], if the Melnikov potential, given by 

00 

L{e)= Y,nh\e)M^\e)), (5) 

t=— 00 

has a nondegenerate critical point in W, then the manifolds W'^{za, fs^e) and W'^i^z^, fs,e) 
intersect transversely for e small enough (the same conclusion is valid for W'^{za., fs^^ 
and W^{zi,, fs^t))- It is known that the series for L converges absolutely and uniformly 
to a C^ function on compact subsets of lA. As an example below, we consider the case 

p[e, e') = cos2 Tie. 

Thus the potential L provides sufficient conditions for the transversal destruction of 
the saddle connection, just as in the classic applications of the Melnikov integral. In §0 
we formulate a slightly stronger version of the previous result (Theorem ^). 

The Melnikov potential also has a physical interpretation: the difference between val- 
ues of L at neighboring critical points is the area of the turnstile lobe that they delineate. 
The area of the lobe gives a coordinate independent measure of the separation of the 
perturbed stable and unstable manifolds as well as the flux from one region to another 
(cf. i). 

The sum for the order e approximation for the lobe area is calculated explicitly in §0 
and compared with numerical results in §^. We summarize our results, as 

Theorem 1 (Main Theorem). Let Ss, given by (0j with potential Vs ([§j, be the gen- 
erating function for the integrable Suris map fs- Then 

Sse = Ss + e COS^ 710 



generates a twist map fs^e- For allQ < 5 <1 and e small enough there are two neighboring, 

•i,0) to Zb= (|,i 



primary heteroclinic orbits from Za = (— |, 0) to Zb = (|, 0). The lobe defined by these two 



orbits has an area given by 

AiS,e) = eTiu)+Oie') (6) 



where z/ = (1 — -\M)/(1 + v^) and 



rfz/. . . 

fe=l 



-i+8E^^t5^ = ''J(°.''). m 



where ^4(2;, g) «s t/ie Jacobi theta function. 

Note that the muhiphers of the hyperbohc fixed point of fs are z/ and l/iy, so that z/ is a 
natural parameter to use (see e.g., |TT1). 

The series r(z/) is rather intriguing. It is an analytic function of u on the interval 
{0 < u < 1}, and approaches zero rapidly as u increases. In fact, though (J^) implies that 
r(z/) is strictly positive, it approaches zero exponentially fast: 

T.. X ^ 47r \' / -7r2 



\og{l/u)J ^\\og{l/iy) 

as z^ ^ 1^, or equivalently 6 -^ 0^. However, note that (^ is only valid for fixed 6 as 
e — > 0. If these parameters both approach zero, for example according to some relation 
6 = S{e), then the analysis becomes much more difficult. In particular when 0(e) term in 
(^) is exponentially small in 6, it could easily be dominated by terms that are formally 
higher order in e. 

2 Melnikov potential for twist maps 

In this section we review the derivation of the Melnikov potential for twist maps [^j . We 
begin with a C^ Lagrangian generating function S{6, 6') that satisfies the twist condition 
and has zero net flux. It gives a map of the annulus implicitly through (^). Alternatively, 
let the action W oi a sequence [6] = {9\ 6^*+^, ... ,9^}, be 



W[9] =J2Si9\9 



t=i 



Then, an orbit of the map that begins at ^* and ends at 9^ corresponds to a critical point 
of W under variation with respect to the interior points 0. The corresponding momenta 
are then defined through (|^) as r* = —diS{9^,9*^^). Thus, for example, a point {9a, To) 
is a fixed point of the map if and only if 9a is a critical point of 5(6', 9), and r^ is defined 
through (^. 

The distance between the perturbed stable and unstable manifolds of any twist map 
with a saddle connection can be computed using the variational principle [|l^, [l^. We 



begin with a map / generated by 5*. Suppose that / has two hyperbolic fixed points 
Za = {Ga^fa) and Zb = {9h,rh), and there is a saddle connection defined by the graph 
of a function x(6') on the interval U = {9a, 9b) between these points. A diffeomorphism 
h : U ^U is induced by the restriction of the map to the saddle connection: 

f{9,x{9)) = {h{9),x{h{9))). (8) 



Let P he a. C^ function with zero net flux. Then the function 

s,ie,e') = s{e,e') + ePie,e') 

generates a twist map /^ for small enough e. Since hyperbolic points are nondegenerate 
critical points of the action |]14|, the perturbed map will have nearby hyperbolic fixed 
points for small enough e. A simple case occurs when 6a is a critical point of P{6, 6) as 
well as of S{6, 9) since it is then a critical point of S^{9, 9) as well. Thus the fixed points 
of S^ will have unchanged configurations, but their momenta will be modified according 
to®. 



The action gives useful formulae for the areas of regions for twist maps |T5[ . We will 
use one such relation to obtain the Melnikov potential: a relation between the graph x 
and the action of orbits on the stable and unstable manifolds of a hyperbolic fixed point 



|T^. Let z^ be a point on the unstable manifold of a fixed point Za that is close enough 
to Za so that the segment of W^i^Za) containing z^ is given by a graph (6',x(6')). Let 
(^*, x(^*)), t < 0, be the preorbit of this z^ . Defining the backward action difference as 



— i 
m^^{9')= J2 [s{9\9'+')-S{9a,9a) 



then the unstable manifold is defined by the graph of the function 

A corresponding formula for the forward action of an orbit on the stable manifold yields 
a formula for the graph of an initial segment of the stable manifold, x^^ of ^fo- 



^{9') = 




'~S{9\9'^^)~~S{9,A) 


x'io) = 







(9) 

The difference between these two actions leads to the Melnikov-like formula for the 
transversal intersection of these manifolds. To summarize: 

Theorem 2. Let S be the generating function for a twist map f that has two hyperbolic 
fixed points Za and Zb with a saddle connection given by r = x(^) for 9 ElA = {9a, 9b). Let 
f induce a diffeomorphism h{9) on the connection. Let S^ = S + eP generate the twist 
map /e, such that the perturbation P has the following properties 

a) P{9a,9a) = P{9b,9b)=0 

b) ^L=..m^) = iL=,,m^) = o 

Then for e > small enough: 

i) the perturbed map has two hyperbolic fixed points near Za and Zb; 



ii) the Melnikov potential 

oo 



(10) 



t=—oo 



converges absolutely and uniformly to a C"^ function on U; and 



iii) if L has a nondegenerate critical point on U, then the unstable and stable manifolds 
of the two fixed points intersect transversely. 

We can relate the Melnikov potential to the area of a lobe in the stable and unstable 
manifolds using the action formula of [l^. Suppose p and q are two neighboring hetero- 
clinic orbits, for example as shown in Fig. 0. Then the area lobe is determined by the 
difference in action between the orbits of p and of q [[1^] 

oo 



;ii) 



t=— oo 



The area A is the signed area below the segment of W^ between q and p and above that 




Figure 2: Resonance for the perturbed Suris map with 8 = 0.5 and e = 0.05. A pair of principal 
homoclinic points are labeled p and g, and the exit and incoming sets are labeled 8 and X. 

Since the action is stationary on an orbit and since the actions of p and q are equal 
when there is a saddle connection, it is easy to see that the area of the lobe is determined, 
to lowest order, by the difference between two critical values of the Melnikov potential L 



(cf. |]T2|, |T3l)- To summarize: 



Theorem 3. Let S^, fe,U, h and L be defined as in Theorem |^. Assume that Op and 6q 
are two points in U such that 



a) L'{d,) = L'ie,) = , 

b) L"{9,) > 0, L'\e,) < 0, 

c) L'{6) ^0 foro^ <e <ep. 

Then the heteroclinic points p and q continue to heteroclinic points of f^, and the stable 
and unstable manifolds of f^ enclose a lobe with area 

A{e) = e{L{e,) - L{e,)) + 0{e') . 

3 The Suris Map 

We will use the results of the previous section with f = fs, S = Ss, 5^ = Ss^e, fe = f5,e 
and the particular perturbation P{6, 9') = cos^ it6. We give formulas for the lobe area of 
the perturbed Suris map. To find the saddle connection, we must find the diffeomorphism 
h induced by the dynamics of the Suris map fs on the connection (H). We will consider 
h defined only on the set U = { — | < 9 < ^}. 
We first give some properties of h. 

Lemma 4. Let h be the diffeomorphism ofU given by 

2 _^ _^ [{^ + ^) tan(f ^) + {v-l] 



KiG) = -arctan r -|— ; , (12) 

where z/ = (1 — v5)/(l + v^) and < 5 < I. Then h^ satisfies 

a) /i^ = hyt, for all t E Z. 

b) h,{-9) = -h,-i{9). 

c) —J is an stable fixed point and ^ is a unstable fixed point of hy. 

d) V^{9) = K{9) -29 + h-\9), for 9eU, for ¥& given by (j^. 

Proof. A direct computation proves a) b) and c). This stability properties are implied by 
^'u{~\) — ^ ^^*i h'v{\) ~ u- Since u G (0, 1) this implies that the former is stable and the 
latter is unstable. For d), we will first show that when 9 eU, 

/I ,+i//ix 2 / ±v5cos7r^ \ ,^„, 

9-hf{9) = -aictan[- 7-—— , (13) 

vr \^1 =F V(5sm7r6'y 

To see this, notice on one hand that ii9 eU, then h{9) G U. Therefore —1 < 9—hy{9) < 1. 
On the other hand, 

t^n (-[9 - h,{9)]) = j^—— . (14) 

V ^ ^ 1 — V sm vrt' 



This implies (|TB|) for the upper signs. The substitution 9 \-^ —9 in equation (0), gives 



the lower signs. It is easy to see that ii 9 G U, then ([1^ is positive and therefore 
< 9 - K{9) < 1. We conclude that for all 9 e U, -1 < 29 - K-i{9) - K{9) < 1. 
To finish the proof we take tangent of the second difference and use the sum formula for 
tangent to obtain 



tan (-[20 



K-ii9) - h,i9)]) = - 



25cos'k9 sinvr^ 

— 5 sin^ Ti9 + 5 cos^ vr^ 
5sm{2'n9) 



1 + 5cos{2ti9) ■ 

Since this last expression is — tan(|V^'(6')), this completes the proof. D 

With the help of lemma ^, we can give a description of the intersection of the saddle 
connection between (— |, 0) and (|, 0). This connection is given by the functions x^, and 
is shown as the heavy curves in Fig. ^ 

Lemma 5. Let fs be the twist map generated by Ss- Then (— i, 0) and (|, 0) are hyperbolic 
fixed points for fg and there are two saddle connections between them given by the graphs 
r = x+i9) =9- K{9), and r = x~iO) = 9- K-i{9). 

Proof. Using Lemma ^, the map can be written in the form 

fs{9, r) = {r + K{9) - 9 + Vi (0), r + K{9) -29 + h,-i (9)) . 
Therefore, ii r = x{9) = 9 — hi,-i{9) then 

fs{9,x-m = {K{e).hM-0) = {K{0),x'{K{e)) ■ 

Now since = — | is an stable fixed point for h,^, this graph gives the left going saddle 
connection. In the same way setting r = x^{9) = 9 — h^{9) gives 

fsiO,x^m = {K-.{9),K-.{9) -9) = {K-.{9),x''{K-m) ■ 

Since the map conjugates to hy-i on this graph, this is clearly the right going saddle 
connection. D 

3.1 Proof of the Main Theorem 

We now can sketch the proof of the main theorem. The analysis of the infinite series for 
the Melnikov potential relies on some summation formulae for elliptic integrals. 

Proof. Let P{9, 9') = cos^ 7t9. It is clear that P satisfies the conditions of Theorem |^. 
Let L be the Melnikov potential, (|^). According to Theorem |^ a sufficient condition for 
transversal intersection of the perturbed manifolds is that L has a nondegenerate critical 
point on the interval W = { — 1<0<|}. The graph of L{9) over U is shown in Fig. ^ 

To proceed, we define L{z) = L(| arctan(z)). A little algebra gives L{z) = J2t^-oo (^t{z) 
where 

4z/2*fl - z^Y 
at{z) 



((l-2)2 + iy2t(l + ^)2)2 



2.8855 - - 



2.8854-- 



2.8852 




0.2 



0.4 



Figure 3: Melnikov potential, L for 5 = 1/2. 



Since z = taii(|^) is monotone increasing, we need only show that there is a nondegenerate 
critical point of L. We are going to show that Zq = is a. local maximum of L and 
Zp = {1 — y/P)/{l + y/u) is a local minimum. 



First, we write L as 

oo 

L{z) = ao{z) + ^{at{z) + a_t(z)} 



(15) 



t=i 



Now, since «<( 
gives 



a^t{z), this implies that L is even, and so L'{0) = 0. Differentiation 



t=i 



128 (1 - 4z/2* + zy4t) j,2t 
(1 + z/2*)^ 



-8ro(z/2) . 



This expression defines a function Fq, which can be rewritten by expanding the denomi- 
nator and rearranging the sums as 



Fo(g) = l + 16^^ 



— l)^k^q 






k=0 



This sum can be obtained from Fourier expansions of Jacobi elliptic functions which are 
written in terms of the nome, q{k) = exp (— tt-^), where K{k) is the complete elliptic 



integral, K' = K{k') and k' = 1 — k'^ |T^. Recall that g is a monotone increasing 



diffeomorphism of k on (0,1). A Fourier expansion for elliptic functions [0, formula 
912.01] can be differentiated to show that 



ro(g) 



77 



k'K] =^l{0,q) 



where ■&4(z,q) is a Jacobi theta function |]T8|, formula 1051.01]. Therefore L"{0) < 0, for 



all < z/ < 1. This shows that 6*^ = is a nondegenerate local maximum of L for all 
< z/ < 1. 



We now wish to show that Zp is the local ininimuni. Rewrite L as 



oo 



Uz) = ^{am W + a-t{z)} • (16) 



t=o 



Since at+i{zp + () = a^t{zp — (), L is even around Zp, and so L\zp) = 0. After some 
algebra, the second derivative is 



fl _L iyl + 2t 



>4 



Notice that L"{zp) = ro(i^^) — ro(z/). Once again, using Fourier series for elliptic functions, 
we find 

ro(g^) = k" ^^K ' 

so that 

u'izp) = 0ir^ (kkr, 

and therefore L"{zp) > 0, for all < z/ < 1. This shows that 6p is a nondegenerate local 
minimum of L, for all < z/ < 1. 

Using Theorem ^ we conclude that a lobe of area 

A(5,6) = 6(L(0,)-L(^p))+O(e2) 

is enclosed by the stable and unstable manifolds, where dp and 6q are given above. Finally, 
we use ( p!5| , |l6D to obtain 

V 



L{e,) = m = i + sj^ji 

t= 

oo 

L{ep)=L{zp) = ^Y.Tl 



.- , Z/2*)2 
°° y2t+l 



jj2t+l\2 ■ 
t=0 ^ ' ^ 

Subtracting these two, and expanding the series we obtain 

\fc L, ,/c 



L{e,) - LiOp) = 1 + 8 5^ \;^, = r(z.) 

fc=i 

where r(z/) was defined in (^. Finally, this result can be explicitly written as 

r(g)=0KA:'^ =^tiO,q), (17) 

upon differentiating the Fourier series [jl8|, 908.51] once. We have therefore obtained 
(I)- □ 

The explicit formula for the first order approximation of the area is compared with 
numerical computations in Sec. 0. 



3.2 Exponentially small behavior 

In this subsection we investigate the asymptotics of the sum, (|^), as 5 — > 0. Our result is 
summarized as: 



Lemma 6. Let r(z/) be defined by (Qj, th 



en 



Air V [ -TT^ 



Proof. We use the expression (|1^) in terms of the eUiptic integral K and the nome q. 



Recall that g ^ 1 as A; — * 1 and k' -^ 0^. The needed asymptotic forms are |]l| 
formulae 112.01 and 112.04] 

K{k) ~ log(i), 



fc' ~ 4 exp 



TT^ 



21og(l/g), 
as k -^ 1^. Putting these into (O) gives the promised result. D 



We again note that ([T8|) does not imply that the area itself is necessarily exponentially 
small: (|^) is only valid for fixed 5 as e — > 0. Nevertheless we will see in §^ that (|1^) agrees 
remarkably well with our numerical calculations of the lobe area. 

3.3 Anti-Integrable Limit 

It is much easier to obtain an expression for the lobe area for large large e. Our main 
purpose is to have an expression that is valid in the regime of large e where numerical 
calculations are difficult. This expression is easy to obtain using the "anti-integrable 



limit" |T9|. For simplicity, in this section we assume that P depends only upon 6 in this 
section. In this case the anti-integrable limit is obtained by scaling the action by e^^, to 
get 

s = sie,e')/e = p{e) + e-^Ss, 

and then setting e~^ = 0. The point is that when e ^ 1 > 6 the points on the two 
heteroclinic orbits are all found in a neighborhood of the critical points of the potential 
P = cos^(7r^), i.e., at 6' = m/2. In the anti-integrable limit, an orbit consists solely of a 
sequence of configuration points sitting at these critical points. Thus the two heteroclinic 
orbits are given by the sequences 

The action is stationary on an orbit, thus to first order in the small parameter e~^, the 
change in the orbit with e can be ignored in the action difference to give 



^ 00 



For the Suris map, this yields 



1 1 

4~^ 



A{5, e) = e - - - — [dilog(l + 5)- dilog(l - 5)] + 0(e-^) 



(19) 



We compare this result with the numerical calculations in §|[ 



4 Numerical comparison 

In this section, we compare the theoretical results with numerical computations of the 
lobe area. The task is to find the actions of the two homoclinic points p and q. For this 
task we use the symmetry of the Suris map. The action difference between the orbits of 
q and p gives the lobe area, (|1T|). 

The Suris map is reversible: it is conjugate to its inverse by an involution R: E? = I 
and RfsR = fs^^, where the reversor is 



R: {e,r) -^ {-e,r +v'{e)) . 



(20) 



Fixed sets, Fix(i?) = {{6,r) : 6 = 0} and Fix{fsR) = {{d,r) : r = 26} of the reversor are 
important because points q e W^{za) nFix(i?) or g G W^{za) nFix{fsR) are heteroclinic 
from Za to Zb [Q. Thus to find a heteroclinic orbit it is sufficient to do a one dimensional 
search for points on the unstable manifold that intersect one of these fixed sets. 



Fi.x(R) 



Fix(fR) 




Fix(f^R) 



Figure 4: Finding some heteroclinic orbits using symmetry. Shown are the fixed sets Fix(/*i?), 
for i = 0, 1, 2 and a rough example of initial guess zq, leading to tc = 1- 

To find these orbits numerically we move along W^{za) to the first points that intersect 
Fix(i?) and Fix(/i?), respectively. For example, to find q G Fixi?, let 2;° = (^°,r°) = 
Za + 77-E" where E^ is the unstable eigenvector, and 77 is a small parameter to be chosen 
below. Let tc + 1 be the first time for which the iterate of z^ is beyond Fix(i?). Now 
choose a point z{6) on the line from z^ to z^ = fs{z^), parameterized by the initial angle 
9. By construction we are guaranteed that the function Z{9) = T^ifl^i^i^^)) has a zero 
for 6 G [6^, 6^]. We use a root finding method (Brent's method) to determine this zero to 
some precision, say p. The choice of precision infiuences the original value for 77, as well 
as the number of iterates until a crossing. Assuming W^ is smooth, the point z^ will be 



0{rf') away from W^. After tc iterates, however, this error will decrease by the factor A"**" 
where A is roughly the unstable multiplier of the fixed point. There is no sense in having 
this error smaller or larger than the precision of our root finder, so we set p ~ X'^^rf. 
On the other hand, since we start a distance t] from the fixed point, and wish to go a 
distance 0(1) to find the first crossing of the symmetry line, we have rjX^" = 0(1). Thus, 
it is appropriate to set i] ~ p^^^. To find the second homoclinic point, p G Fix{fsR), we 
repeat the above analysis, using crossing of Fix{fsR) to determine tc, etc. The lobe area 
is given by the difference in action between these two orbits, from (pT]). 

For our computations, using IEEE double precision arithmetic, we set p = 10^^^. 
These computations give apparently accurate results providing A » 10~^^. Subsequent 
to our obtaining these results, Delshams and Ramirez-Ros used extended precision arith- 
metic to obtain lobe areas as small as 10"*^^°° f|T|. 
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Figure 5: Log-Linear plot of Aje as a function of b for various values of e compared with the 
theoretical expression (0) (thick line) and the small b expression (O) (thin line) 



In Fig. 1^ we show a comparison of the result of (|^) with the numerical results on a 
log scale. The analytical result agrees well with the numerical results when e = 10"^. We 
show the same data on a linear scale in Fig. ^. We plot the sum (^ only for 5 > 0.15, 
since it is indistinguishable from asymptotic formula (|l^) for smaller values. Moreover, 
the asymptotic formula agrees with the e = 10~^ computation within 1% up to 5 = 0.8. 
In Fig. ^ the anti-integrable results evaluated at e = 1 are also shown. We are unable 
to obtain numerical results for such a large e, as the multiplier of the fixed points is too 
large. 
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Figure 6: Linear plot of Aje as a function of b for various values of e. The solid curve is the 
anti-integrable limit (19) for e = 1. 



5 Conclusion 

The perturbed Suris map studied here depends on two parameters, the Suris parameter 5 
and the perturbation strength e. We obtained the lobe area for the fixed point resonance 
of this map for small e to 0(e) and for large e to 0(e~^). In the small epsilon case, the 
numerical computations indicate that the asymptotic formula (|18D for 5 -^ is a good 
approximation even up to 5 = 0.8. 

A similar Melnikov analysis is possible for other standard maps [0, |lT|. In particular 
01 showed that there exists a large class of standard maps with saddle connections. In- 
terestingly these maps are not all integrable. The use of the Melnikov potential is not 
restricted to twist maps; in fact, it can be applied to any exact symplectic map [p!3[ . 

The method can also be applied to any higher dimensional twist map that has a saddle 
connection of the type described in this paper 



see 



13| , (T^l)- Integrable examples of such 



maps have been found in |^, ^ ^, and the Melnikov method has been applied to applied 
to the four dimensional McLachlan map [1^. Similar techniques are also applicable to 
volume preserving maps [24|. The study of perturbations of twist maps with saddle 



connections in higher dimensions is important because it could help in the development 
of a higher dimensional theory of transport. 



References 

[1] A. Delshams and R. Ramirez-Ros. Singular separatrix splitting and the Melnikov 
method: an experimental study. Experiment. Math., 8(l):29-48, 1999. 

[2] Yu. B. Suris. Integrable mappings of the standard type. Func. Anal, and AppL, 
23:74-76, 1989. 



[3] J.D. Meiss. Symplectic maps, variational principles, and transport. Rev. Mod. Phys., 
64(3):795-848, 1992. 

[4] V.K. Melnikov. On the stability of the center for time periodic perturbations. Trans. 
Moscow Math. Soc, 12:1-57, 1963. 

[5] R.W. Easton. Computing the dependence on a parameter of a family of unstable 
manifolds: Generalized Melnikov formulas. Non. Anal. Th. Meth. and AppL, 8:1-4, 
1984. 

[6] M.L. Glasser, V.G. Papageorgiou, and T.C. Bountis. Melnikov's function for two- 
dimensional mappings. SIAM J. Appl. Math, 49:692-703, 1989. 

[7] A. Delshams and R. Ramirez-Ros. Poincare-Melnikov-Arnold method for analytic 
planar maps. Nonlinearity, 9(1): 1-26, 1996. 

[8] H. Lomeli. Perturbations of elliptic billiards. Phys. D, 99:59-80, 1996. 

[9] H. Lomeli. Saddle connection and heteroclinic orbits for standard maps. Nonlinearity, 
9:649-668, 1996. 

[10] H. Lomeli. Applications of the Melnikov method to twist maps in higher dimensions 
using the variational approach. Erg. Th. and Dyn. Sys., 17:445-462, 1997. 

[11] A. Delshams, R. Rami rez Ros, and T.M. Seara. Splitting of separatrices in Hamil- 
tonian systms and symplectic maps. In C. Simo, editor, Hamiltonian Systems with 
Three or More Degrees of Freedom, volume 533 of Series C: Mathematical and Phyis- 
cal Sciences, pages 39-54. Kluwer Academic, 1999. 

[12] R.S. MacKay and J.D. Meiss. Relationship between quantum and classical thresholds 
for multiphoton ionization of excited atoms. Phys. Rev. A, 37(12):4702-4706, 1988. 

[13] A. Delshams and R. Ramirez-Ros. Melnikov potential for exact symplectic maps. 
Comm. Math. Phys., 190(l):213-245, 1997. 

[14] R.S. MacKay and J.D. Meiss. Linear stability of periodic orbits in Lagrangian sys- 
tems. Phys. Lett. A, 98:92-94, 1983. Relation between Residue and second variation 
of action. 

[15] R.S. MacKay, J.D. Meiss, and I.C. Percival. Resonances in area preserving maps. 
Phys. D, 27:1-20, 1987. 

[16] E. Tabacman. Variational computation of homoclinic orbits for twist maps. Phys. 
D, 85:548-562, 1995. 

[17] R.S. MacKay, J.D. Meiss, and I.C. Percival. Transport in Hamiltonian systems. Phys. 
D, 13:55-81, 1984. 

[18] P.P. Byrd and M.D. Friedman. Handbook of Elliptic Integrals. Springer- Verlag, 1971. 



[19] R.S. MacKay and J.D. Meiss. Cantori for symplectic maps near the anti-integrable 
limit. NonUneanty, 5:149-160, 1992. 

[20] D. Sterling, H.R. Dullin, and J.D. Meiss. Homoclinic bifurcations for the Henon map. 
Physica D, 134(2): 153-184, 1999. 

[21] R.I. McLachlan. Integrable four- dimensional symplectic maps of standard type. Phys. 
Lett. A, 177(3):211-214, 1993. 

[22] J. Moser and A. Veselov. Discrete versions of some classical integrable systems and 
factorization of matrix polynomials. Commun. Math. Phys., 139:217-243, 1991. 

[23] Yu. B. Suris. A discrete-time Gamier system. Phys. Lett. A, 189:281-289, 1994. 

[24] H. E. Lomeli and J.D. Meiss. Heteroclinic primary intersections and codimension 
one Melnikov method for volume preserving maps. Chaos, in press, 1999. 



